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ABSTRACT 


The  salt-gradient  solar  pond  was  modelled  both  experi¬ 
mentally  and  numerically.  The  experimental  model  used  a  small- 
scale  tank  with  artificial  sunlight.  Both  temperature  and  con¬ 
centration  measurements  were  taken.  The  numerical  model  used 
a  one  dimensional  heat  conduction  model  to  describe  the  heating 
phenomena  in  the  pond. 

The  two  models  predicted  the  same  type  of  behavior 
as  that  expected  in  full-scale  ponds.  There  were  small  dif¬ 
ferences  in  the  two  models.  These  differences  arose  from 
simplifications  in  th6  numerical  model. 
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1.  INTRODUCTION 
1 . 1  Definition 

The  salt-gradient  solar  pond  is  a  body  of  water  that 
is  used  as  both  collector  and  storage  of  solar  energy.  Salt  is 
added  to  the  pond  such  that  the  concentration  of  salt  increases 
with  depth,  thus  inhibiting  free  convection.  The  pond  can  be  from 
2i  to  10  meters  deep,  with  fresh  water  at  the  surface  and  a 
nearly  saturated  solution  at  the  bottom.  The  bottom  of  such 
a  pond  is  usually  black,  for  high  absorption.  When  the  water 
in  the  bottom  region  of  the  salt-gradient  pond  is  heated,  it 
tends  to  remain  hotter  than  the  rest  of  the  pond,  as  convection 
is  suppressed.  Thus,  heat  extraction  is  made  from  the  bottom 
of  the  pond. 


1 . 2  Previous  investigation 

The  first  serious  research  effort  in  the  area  of 
artificially  created  solar  ponds  was  begun  in  1958  by  H.  Tabor 
(Ref.  1).  This  work  continued  until  1966  when  the  interest  in 
solar  energy  decreased  due  to  the  low  cost  of  fuel -oil.  After 
the  energy  crisis  of  1973,  interest  in  solar  ponds  awakened. 

The  Scientific  Research  Foundation  in  Israel  began  solar  pond 
research  in  1974.  Two  demonstration  power  plants  have  been  set 
up  in  Israel,  both  of  which  operate  at  temperatures  near  90°C. 
Similar  research  was  begun  in  1974  at  Ohio  State  University  in 
the  U.S.  Most  of  the  experimental  work  reported  from  both 
Israel  and  the  U.S,  has  been  performed  on  full-scale  ponds, 
with  surface  areas  ranging  from  150  ra^  to  7000  m^. 

Some  laboratory  experiments  have  been  reported, 
however,  which  investigate  the  concentration  and  temperature 
gradients,  using  tanks  of  0.5  m  diameter  (Ref.  2).  These  tanks 
were  heated  from  the  bottom  and  thus  did  not  simulate  solar 
radiation. 
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Wilkins  and  Pinder  have  reported  work  done  on  a  small- 
scale  solar  pond  model  (Ref.  3).  The  pond  they  used  had  surface 
dimensions  of  65  cm  65  cm  and  was  25  cm  deep.  Their  model 
used  partitions  to  separate  the  non-convecti ng  layer  from  the 
convecting  layers,  and  was  exposed  to  actual  radiation  in 
Vancouver. 


In  1981,  investigators  at  VKI  attempted  to  model  the 
"saturated"  solar  pond,  with  radiation  supplied  by  five  halogen 
lamps  (Ref.  4). 

The  present  study  differs  from  the  others  in  that  a 
non-convecti ng  solar  pond,  without  partitions,  was  modelled  on 
a  small  scale  with  controlled  artificial  sunlight.  The  experi¬ 
mental  results  are  then  compared  to  the  results  from  a  numerical 
model . 


1 .3  The  motivation  for  laboratory  modelling 

If  solar  ponds  are  to  be  used  for  energy  production, 
a  model  is  needed  to  test  various  designs  and  improvements. 

A  typical  solar  pond  has  a  surface  area  on  the  order  of  1000  m^. 
Thus  it  would  be  highly  impractical  to  build  a  full-size  solar 
pond  to  test  any  new  idea  or  possible  improvement.  There  are 
many  existing  computer  simulations  of  solar  ponds 
(Refs.  5, 6. 7, 8). 

Numerical  models  are  limited,  however,  to  the  inclu¬ 
sion  of  phenomena  that  can  be  described  mathematically.  Although 
it  may  be  argued  that  this  is  always  possible,  it  is  certainly 
not  always  practical.  Thus  the  need  for  reliable  small-scale 
laboratory  models  becomes  evident. 
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2.  THEORY 

2 . 1  Basic  principle 

Solar  radiation  is  received  by  the  pond  through  the 
pond  surface.  Some  of  the  incoming  radiation  is  reflected  by 
the  surface  and  some  is  absorbed  near  the  surface  (see  Fig.  1). 
Much  of  the  radiation  reaches  the  bottom  of  a  shallow  (<  10  m) 
pond  and  is  absorbed  by  the  absorbing  bottom.  Figure  '  hows  the 
amount  of  transmitted  radiation  as  a  function  pond  d  h.  It 
can  be  seen  from  this  figure  that  40%  of  the  incoming  ^oiation 
can  reach  the  bottom  of  a  1.5  m  pond.  Thus,  the  bott  "uld 
absorb  the  heat  and  consequently  re-radiate  it,  heatii  ,.ie  pond 
from  the  bottom.  In  a  pure  water  pond,  the  hot  water  would  be 
lighter  than  the  colder  surface  water  and  free  convection  would 
begin,  mixing  the  water  and  allowing  the  heat  to  escape  into 
the  atmosphere  through  the  surface.  In  a  salt-gradient  pond, 
however,  there  exist  dissolved  salts  in  the  pond  in  proportions 
that  increase  with  pond  depth.  That  is,  the  solution  at  the 
pond  surface  contains  very  little  salt  and  that  at  the  pond 
bottom  is  nearly  saturated.  This  increase  in  solution  density 
conteracts  the  decrease  in  density  due  to  the  temperature  change, 
thus  keeping  the  hot  solution  at  the  bottom  of  the  tank.  The 
stability  requirement  can  be  expressed  as  : 

ifi.  =  fie.)  i£.  +  fie.]  ia  >  0  (i) 


3Tjq  ah 


aqjj  ah 


where  : 


density  of  pond  solution  [=1  kg/m^ 

distance  below  the  surface  (positive  downwards)  [ »1  m 
temperature  I  *1  “C 

salt  concentration  ( *1  ^ 

-3 


I  t 


i 
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As  long  as  the  criteria  expressed  in  equation  (1)  is  maintained, 
the  pond  will  be  non-convecting  and  the  heating  phenomena  will 
be  heat  conduction. 


2 . 2  The  convective  layers 


If  the  entire  pond  were  truly  non-convecti ve ,  the 
temperature  gradient  would  be  as  shown  in  figure  3a.  The 
gradients  circled,  however,  quickly  become  great  enough  to 
overcome  the  suppressive  action  of  the  density  gradient  and  con¬ 
vection  begins  in  the  zones  shown  in  figure  3b.  Thus  in 
practice,  a  salt-gradient  pond  is  composed  of  three  layers  : 

(1)  The  upper  convective  zone  (UCZ). 

(2)  The  non-convecti ve  zone  (NCZ). 

(3)  The  lower  convective  zone  (LCZ), 

The  parallel  temperature  and  density  gradients  occur  In  the  non- 
convecting  zone,  as  shown  in  figure  4. 

2.3  The  non-dimensional  equations  (Ref.  5) 


A  heat  balance  over  a  small  layer  in  the  non-convecting 
zone  yields  the  following  equation  : 


il 

at 


a 


ill' 

3h2 


+ 


sec 


(2) 


where  : 


t  =  time  [  =]  s 

k  n  ^ 

a  =  thermal  diffusivity  =  -  [  =]  — 

oCp  s 

w  =  effective  attenuation  coefficient  =  m 

C  =  heat  capacity  of  solution  [ =]  — — 

P  kg»C 
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I =  radiation  at  a  given  pond  depth 

=  I  (l-F)exp{-u[  h-61  sece  ]  t =1  — 


0 

r 

n 

f 

5 

I 


refracted  angle  =  sin 


sine^ 

n 


I  1  rad 


solar  zenith  angle  [  =]  rad 

refractive  index  of  solution  [=1  non-dimensional 

fraction  of  radiation  absorbed  in  a  small  distance,  6, 
of  solution,  non  dimensional 

small  distance  in  which  fraction  F  is  absorbed  t =]  m 


irradiance  just  beneath  the  pond  surface  [ =] 


In  equation  (2)  it  is  assumed  that  heat  transfer  occurs  in  one 
direction  only.  The  depth  of  a  solar  pond  is  ordinarily  much, 
much  smaller  than  either  the  length  or  the  width,  making  this 
assumption  good  for  such  an  application.  Also,  it  is  assumed 
that  all  the  radiation  reaching  the  bottom  is  absorbed  by  the 
bottom.  An  equivalent  situation  occurs  when  some  radiation  is 
reflected  from  the  bottom,  in  diffuse  form,  but  is  consequently 
absorbed  in  the  storage  zone  (the  lower  convecting  zone). 

The  solution  of  eq.(2)  requires  one  initial  condition 
and  two  boundary  conditions.  The  initial  condition  can  be  chosen. 
One  boundary  condition  evolves  from  a  heat  balance  over  the  upper 
convecting  zone 


at 


phiC, 


|^l-(l-F)exp{  -  u(h-6  )sec6j,}J 


+ 


_k _ 

phjCp 


faT 


ah 


h=hi 


(3) 


where  : 
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=  temperature  in  the  upper  convecting  zone  [ =]  "C 

hj  =  depth  of  the  upper  convecting  zone  [  =]  m 

=  rate  of  energy  loss  at  the  surface  [ =]  W/s 

The  first  term  on  the  right  represents  the  heat  absorption  in 
the  UCZ  due  to  radiation.  The  second  term  represents  the  heat 
conduction  across  the  boundary  between  the  NCZ  and  the  UCZ. 

The  last  term  represents  heat  lost  through  the  pond  surface. 

The  boundary  condition  at  the  bottom  can  be  found  by  performing 
a  similar  heat  balance  over  the  LCZ. 


— ^  ^  ^  I  exp  I -u(D-6  )  sece 

at  pD^Cp  s  L 


k  far' 

pOgCp 


pD^C 


p 


Q 

pOgC 


(4) 

(4) 


where  : 


Tjj  =  temperature  in  the  LCZ  (  =]  "C 
0^  =  depth  of  the  LCZ  [ =1  m 

D  =  hi  +  depth  of  the  NCZ  [ =J  m 
k  =  thermal  conductivity  of  sorn[=]  W/m“C 
Qlg  =  rate  of  energy  lost  to  the  ground  [=1  J/m^ 

Qq  =  rate  of  energy  transfer  to  load  [  =]  J/m^ 

The  first  term  on  the  right  represents  the  heat  absorption  in 
the  LCZ.  The  second  term  represents  the  heat  conduction  across 
the  boundary  between  the  LCZ  and  the  NCZ.  The  third  term 
represents  the  losses  to  the  ground  and  the  last  term  represents 
the  heat  withdrawn  from  the  pond  to  meet  the  load. 
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To  non-dimensional ize  equations  2,  3  and  4,  it  is 
necessary  to  choose  a  characteristic  time,  temperature  and 
distance.  The  characteristic  temperature  is  chosen  as  simply 
Tq,  a  refe'^v'nce  temperature,  which  can  be  the  initial  pond 
temperature,  or  the  ground  temperature,  or  in  any  other  con¬ 
venient  way.  The  characteristic  distance  is  chosen  as  D, 
the  combined  depth  of  NCZ  and  UCZ.  The  characteristic  time 
is  chosen  as  which  can  be  thought  of  as  a  "thermal 
travelling  time"  over  the  distance  D.  Thus  equations  2,  3  and 
4  become,  in  dimensionless  form. 


3^  _  3^  ^  D^I^tisece^ 
3t  3Z2  IcTo 


3t  kToh^ 

.  D 


j^l-(l-F)exp  { -w(hi-5)sece^}J 


3e 


39 


b  _ 


3t 


hi  i3Z 

DM1-F)l3 

kToD, 


Z  =Zi 


khi 


1  -  - 

f  \ 

39 

J  D 

e 

,3Z. 

Z  =  1 


kD 


("b-'Gw) 


V! 

kToD 


where  : 


T  .21 

D2 


(5) 


(6) 


(7) 


Z 


h 
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0  = 

S  T 

T  0 


U|_Q  =  heat  transfer  coefficient  at  pond  bottom  [  =]  W/m^°C 
=  heat  transfer  coefficient  at  surface  [  =]  W/m^^C 

Q  ^ambient 

®a - - - 

1  0 

Aground 
T. 


'GW 


Several  dimensionless  "numbers",  or  parameters, 
appear  in  the  formulation  of  the  above  equations  : 

f 

(1)  The  Fourier  number  =  Fq  =  —  where  t  =  characteristic 

D2 

^  n2 

time  (  =]  LC  (t  ^ —  here,  in  order  to  choose  a  t  to  maintain 

a 

the  physical  meaning  of  the  numbers).  Fq  gives  the  magnitude 
of  the  contribution  to  the  change  in  temperature  in  the  NCZ 
by  conduction. 


(2)  N1  = 


(3)  N2  = 


(4)  N3 


aCpTo 


_ 

PhiCpTo 


the  magnitude  of  the  contribution  to  the 
change  in  the  NCZ  by  radiation 


the  magnitude  of  the  contribution  to  the 
change  in  temperature  in  the  UCZ  due  to 
radiation 


kt 


-  the  magnitude  of  the  contribution  to  the 
Phi  C  D  change  in  temperature  in  the  UCZ  due  to 

P  conduction 


tu 


(5)  N4 


LS 


phiC, 


the  magnitude  of  the  contribution  to  the 
change  in  temperature  in  the  UCZ  due  to 
losses  through  the  surface 


I 
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(6)  N5  =  -  =  the  magnitude  of  the  contribution  to  the 

TqpD  C  change  in  temperature  in  the  LCZ  due  to 

^  P  radiation 


(7)  N6  =  -  =  the  magnitude  of  the  contribution  to  the 

pD  C  D  change  in  temperature  in  the  LCZ  due  to 
^  P  conduction 


Uigt 

(8)  N7  =  — =  the  magnitude  of  the  contribution  to  the  change 
pD  C  in  temperature  in  the  LCZ  due  to  losses  to 
^  P  the  ground 

In  this  formulation,  we  will  not  consider  the  load  term,  as  this 
study  is  concerned  with  the  behavior  of  a  solar  pond  with  no 
heat  load.  For  a  typical  solar  pond  of  total  depth  1.5 
meters,  the  following  values  can  be  calculated  for  the 
dimensionless  numbers  : 


Fo  10 
N1  a  10‘ 


N2  a  10' 
N3  as  lO' 


N4  as  10 


N5  as  1 
N6  as  io‘ 


N7  as  10 


Since,  in  the  NCZ,  Fq  as  N1 ,  both  the  magnitude  of  N1  and  Fg  must 
be  respected.  Since  in  the  UCZ,  N2  »  N3  and  N2  »  N4,  only  N2 
need  be  modelled.  Similarly,  in  the  LCZ,  N5  »  N6  and  N5  »  N7, 
so  N6  and  N7  can  be  sacrificed.  Thus,  four  numbers  must  be 
respected  in  any  experimental  model  :  Fg,  N1 ,  N2  and  N5. 


J 


»  a  *4.  A  *.-! 
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2.4  Oi scretization  of  the  equations 


The  non-dimensional  equations  5-7  must  be  discretized 
for  solution  on  the  VAX  computer.  A  Crank-Nichol sen  scheme 
(implicit/explicit)  was  used.  Using  a  forward  time  discretiza¬ 
tion  and  centered  space  discretization  equation  5  becomes  : 


n+1  n 


-®i 


At 


n+1 


n+1  n+1 


®i+r29i 


"®i-i 


L 

AZ2  J 

r  n  n  n  -| 

1  1 

r ,  n 

1 

I 

D2l^usece^| 

2 

AZ2 

^  i 

- 1 

o 

1- 

_ 1 

(8) 


Using  a  forward  discretization  in  both  time  and  space,  equa¬ 
tion  6  becomes  : 


n+1  n 

~  n+1 

n+r 

•Vi 

n  n 

®i  _  1  0 

02 

-01 

+  1  D_ 

02-01 

At  2  hi 

L  AZ  _ 

2  hi 

_  aZ  . 

D^l 


kToh; 


l-(l-F)exp{-u(hi-6)sec0^} 


“ls 


kh 


1 

- 01 

2  kh 


2  kh  ^ 


(9) 


Using  a  forward  time  discretization  and  a  backward  space 
discretization,  equation  7  becomes  : 
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f  n+1 

n 

■  n+1  n+1' 

n  n  ' 

1 
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j  At 
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aZ 

■  2  D3 
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+ 
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(Note  that  the  heat  withdrawal  term  has  been  dropped) 

Equations  8,  9  and  10  can  be  programmed  in  Fortran 
to  set  up  a  tridiagonal  matrix  to  solve  for  e  in  a  time-marching 
style  (see  Appendix  1  for  the  program  listing).  The  Crank- 
Nicholsen  scheme  is  unconditionally  stable. 

2 . 5  Temperature  behavior 

As  the  pond  is  heated  at  a  constant  heat  input  one  would 
expect  an  exponential  increase  in  temperature.  When  cycling  occurs, 
the  average  heat  input  (average  with  respect  to  time)  is  decreased 
and  thus  the  maximum  temperature  would  be  less  that  that 
with  the  same  constant  heat.  Unequal  cycling,  with  the  time  'on' 
greater  than  the  time  'off  would  result  in  an  intermediate  maximum 
temperature . 


i 

i 


f 
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3.  EXPERIMENT 

3 . 1  Apparatus 

3.1.1  The  £orid 

The  model  used  for  these  experiments  had  a  square 
surface  area  of  0.8231  m^.  The  metal  tank  was  lined  on  the 
inside  wall  with  4  cm  of  glass  bead  insulation  (see  Appendix  2). 

A  pond  depth  of  30  cm  was  used  with  a  mylar  sheet  floating  on  the 
surface.  The  purpose  of  the  surface  sheet  was  to  prevent 
excessive  losses  due  to  evaporation.  A  photograph  of  the  pond 
appears  in  figure  5. 

3.1.2  T  h£  £uji 

Four  halogen  lamps  were  used  to  supply  artificial 
sunlight  to  the  pond.  They  were  situated  at  a  height  that 
provided  energy  to  the  pond  somewhat  uniformly  and  of  magnitude 
(at  pond  surface)  of  550  W/m^.  Figure  6  shows  the  radiation 
Isolines  at  the  pond  surface.  The  standard  deviation  of  the 
data  was  45.7  W/m^,  which  is  8.2%  of  the  average.  For  the 
cycling  phase  the  "sun"  was  controlled  by  a  microprocessor 
clock. 


3.1.3  Tern p£r£t^r^  me£S£reroen_t  s 

The  temperature  of  the  pond  was  measured  by  twelve 
Copper-Constantan  thermocouples  at  the  positions  shown  in  Table  1. 
The  temperature  measurement  chain  is  shown  schematically  in 
figure  7.  A  mechanical  scanning  device,  the  "clop-clop",  was 
implemented  to  scan  the  thermocouples  every  14  minutes.  An 
electronic  T-reference  was  used  that  also  provided  amplification 
for  the  signal  to  be  received  by  the  data  logger.  The  data 
logger  recorded  the  amplified  signal  on  magnetic  tape  to  be  later 
processed  by  the  POP  computer. 
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3.1.4  Salt  concentration  measurements 

A  probe,  operating  on  optical  principles,  was  developed 
to  measure  the  concentration  of  MgCA2  in  the  pond.  A  schematic 
of  the  measurement  chain  is  shown  in  figure  8.  A  laser  beam  is 
directed  into  a  1  mm  plastic  optical  fiber  by  a  focusing  device 
(see  Fig.  9  for  detail).  The  optical  fiber  then  makes  a 
U-turn.  This  U-turn  is  exposed  to  the  solution  in  which  the 
concentration  is  to  be  measured.  As  the  laser  beam  hits  the 
fiber  wall,  some  of  the  light  is  refracted  into  the  solution 
and  some  is  reflected  back  into  the  fiber.  The  amount  of  light 
that  is  refracted  into  the  solution  is  a  function  of  the  index 
of  refraction  of  the  solution  which,  in  turn,  is  a  function  of 
salt  concentration.  Thus,  a  measure  of  the  light  intensity 
change  through  the  bend  is  a  measure  of  the  salt  concentration. 

To  measure  the  intensity  of  light  leaving  the  solution,  a 
photo-resistor  is  used  in  the  receiving  network  as  shown  in 
figure  10,  and  the  output  is  read  in  volts.  Seven  probe  shapes 
were  tested  (Fig.  11).  The  final  shape  chosen  was  probe 
number  7,  whose  calibration  curve  is  shown  in  figure  12. 


3 . 2  The  procedure 


3.2.1  Fil^l2n£ 


The  pond  was  filled  using  a  port  at  the  bottom  of  the 
tank.  It  was  filled  in  layers,  adding  40  i  fresh  water  first. 

gMgCij 

Then  40  i  of  0.05  -  solution  was  added,  very  slowly,  so 

mt  soln 

I) 

as  not  to  mix  the  two  layers.  The  concentrations  were  increased, 
as  shown  in  figure  13,  until  the  last  layer  of 
gMgCti 

0.3  -  was  added.  The  last  layer  was  12  cm  deep  to  induce 

mt  soln 

an  effective  LCZ.  During  filling,  a  concentration  scan  was 
performed  to  monitor  the  filling  method. 
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3.2.2  iie£t2n£ 

The  lights  were  turned  on  and  the  pond  was  heated  for 
75  hours.  Data  was  taken  intermittently  during  this  period. 
Concentration  scans  were  performed  during  the  test  and  at  its 
compl etion . 

3.2.3  Cool 2n£ 

After  75  hours  of  heating,  the  lights  were  shut  off 
and  the  pond  was  allowed  to  cool  for  70  hours.  During  the 
cooling  phase  temperature  data  was  taken  automatically  every 
14  minutes.  At  the  end  of  the  test,  a  concentration  scan 
was  made. 


3.2.4  Da^/ni3.ht  £ycl2n£ 

To  study  the  cycling  in  the  steady-state  mode',  the  pond 
was  first  heated  for  several  days  until  a  surface  temperature  , 
of  about  40*C  was  reached.  Then  the  controlling  clock  was  set 
for  4.8  hours  on,  4.8  hours  off.  A  run  of  16  hours  was  taken. 
Then  the  clock  was  set  for  4.8  hours  on,  2.93  hours  off,  and 
60  hours  of  data  were  recorded.  During  all  cycling  phases, 
temperature  data  was  taken  automatically  every  14  minutes  on 
magnetic  tape.  At  the  completion  of  the  cycling  tests,  a 
concentration  scan  was  made. 

3 . 3  Uncertainty  analysis 

3.3.1  Teinper£t£r£  measurements 

3. 3. 1.1  Est1raation_of  uncertainty 

The  temperatures  were  computed  using  the  following 

formula  : 

T  .  VB2-4AC 


ZA 
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where  : 

B  =  3.03  ±  0.01 
A  =  0.00641  ±  0.00001 
C  =  -  DIG  ±  10 

(where  OIG  =  the  numerical  value  recorded  on  tape). 

Thus,  the  uncertainty  can  be  estimated  using  the  following 
formula  : 


6T 


2 

2  / 

2 

2 

6B 

+ 

+ 

fDIG 

13B  J 

[3A 

Uc  J 

S  DIG 

For  T  =  60“,  DIG  =  500  and 


5T  =  2“ 


Thus,  the  uncertainty  is  6%. 

3. 3. 1.2  Radiation  effects 

The  most  important  thermocouple  error  would  most 
likely  be  from  radiation.  Figure  14  compares  two  sets  of  data 
taken  at  the  same  time.  One  was  taken  with  the  lights  off  and 
the  other  with  the  lights  on  to  observe  the  radiation  effects 
on  the  thermocouple  measurements.  It  can  be  clearly  seen  that 
the  effects  of  radiation  are  felt  only  at  the  first  thermocouple 
near  the  surface.  From  these  data  it  was  decided  not  to  include 
a  correction  for  radiation  effects. 


t 


\ 


I 


3.3.2  £o]ic en t r a t i^o n  inea si^r^me n t s 

The  final  probe  tested  has  a  repeatability  of 
±  0.015  volts.  The  calibration  procedure  has  a  great  influence 
on  the  uncertainty.  A  linear  formula  is  used  for  the  tempera¬ 
ture  variation  with  volts. 

C  =  mV  +  b 

If  the  standard  solutions  are  accurate  to  ±  5%,  one  can  assume 
that  m  and  b  are  accurate  to  about  5%  also. 

m  =  -  0.212  ±  0.010 
b  =  0.239  ±  0.012 

The  uncertainty  in  C  can  then  be  calculated  using  : 


Using  this  formula  for 

C  =  0.22 
V  =  4.16 


one  obtains  : 


17 


4.  THE  NUMERICAL  MODEL 

Equations  3,  9  and  10  (presented  in  section  2.4) 
were  programmed,  in  FORTRAN,  to  be  run  on  the  VAX  computer. 
First  the  equations  were  solved  with  no  radiation  in  the 
program  DECAY  (see  Appendix  I).  Then,  in  the  program  HEAT 
(Appendix  I),  the  radiation  effects  were  added.  Then,  the 
two  programs  were  combined  in  SWITCH  (Appendix  I),  to  enable 
the  day/night  cycling  to  be  simulated. 

Runs  were  made  that  compared  with  the  experimental 
heating,  cooling  and  cycling.  For  the  experimental  comparison, 
the  following  parameters  were  used  : 


lights  "off" 
lights  "on" 


hj  =  0.4  cm 
D  =  18  cm 
Dj  =  12  cm 
u  s  0. 1  m  ^ 
I5  =  550  W/m^ 
At  =  0.00397 
AZ  =0.1 


“ls  '  0 


1I.125  U/m^'c 


LG 


I 
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5.  ANALYSIS  OF  RESULTS* 

5 . 1  Experimental 

5.1.1  Hea  ti^ng 

5 . 1 . 1 . 1  Data 

The  series  of  profiles  shown  in  figure  15  are  selected 
from  the  data  taken  during  the  heating  phase.  Curve  1  repre¬ 
sents  the  initial  conditions  of  the  pond,  that  is,  nearly 
uniform  temperature.  The  uppermost  data  point  of  each  curve 
is  30  cm  (pond  surface)  and  the  next  point  is  at  27  cm 
(see  Table  1).  Thus,  one  can  say  from  data  presented  here,  that 
the  UCZ  is  less  than  3  cm  thick.  That  is, 

0  <  hi  <3  cm 

The  data  shows  nearly  uniform  temperatures  for  the  bottom  12  cm. 
Thus , 

0^  3:  12  cm 

A  pronounced  temperature  gradient  is  formed  during  the  first 
8  hours  of  heating,  representing  exaggerated  heat  absorption 
in  the  first  zone.  In  the  subsequent  curves,  the  gradient 
becomes  less  pronounced,  as  the  heat  is  convected  through  to 
the  lower  regions. 

Figures  15a, b,c  represent  graphs  of  the  data  taken  during 
heating  that  express  the  temperature  as  a  function  of  time  of 
three  of  the  12  thermocouple  positions.  Each  curve  represents 
an  exponential  increase  with  time.  One  can  see  the  time 
constant  is  smaller  near  the  pond  surface. and  larger  near  the 

*  The  complete  set  of  raw  data  is  stored  on  floppy  disk  and  is 
available  on  request. 
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pond  bottom.  In  the  series  of  figures  the  effect  of  radiation 
coming  in  at  the  window  can  be  clearly  seen.  There  are  three 
"humps"  in  the  data  that  can  be  observed  at  t  «  6  hours, 
t  =«  30  hours  and  t  «  55  hours.  These  times  correspond  to  late 
afternoon  from  15  to  17h  on  each  of  the  three  testing  days. 

Figure  17  is  the  same  time  history  for  the  average 
temperature  during  heating.  The  three  "humps"  are  clearly 
visible  here  also. 

5. 1.1. 2  Heat_storage_during  the  heating  phase 


An  analysis  of  the  pond  heat  storage  during  the  heating 
phase  can  be  made  from  the  data  shown  in  figure  17.  The  heat 
input  was  calculated  by  the  known  radiation  at  the  pond  surface. 
The  energy  stored  in  the  pond  was  calculated  by 


heat  stored  =  mCpAT(t} 


where  : 


m  =  mass  of  pond  =  293  kg 

Cp  =  heat  capacity  of  pond  =  3634  J/kg®C 

AT(x)  =  temperature  difference  from  t  =  0  to  t. 


In  such  a  manner,  figure  18  was  generated.  The  minimum  storage 
appears  to  be  about  30%. 

5. 1.1. 3  Global_heat  loss_coef f i cient 

A  global  heat  balance  on  the  pond  yields  : 


mC_  ^  =  I  S  -  HS(T-T  ) 
P  dt  ^  a 


where 


S  =  pond  surface  area  =  [ml 
H  s  global  heat  loss  coefficient  [ =]  W/m^^C 
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Using  equation  (11),  one  can  solve  for  H(T-T  ),  The  variation 

a 

of  HAT  is  shown  in  figure  19.  As  expected,  HAT  increases  with 
time,  that  is,  as  the  average  temperature  increases. 

5.1.2  £ooli^n£ 

Temperature  profiles  taken  during  the  cooling  phase 
are  shown  in  figure  20.  Curve  1  is  the  initial  condition,  the 
heated  tank.  In  the  first  9  hours,  a  drastic  change  takes  place 
as  the  temperature  gradient  changes  sense.  Heat  is  lost  quickly 
at  the  pond  surface.  After  54  hours  of  cooling  (curve  7), 
there  is  still  heat  storage  in  the  LCZ. 

The  time  variations  of  the  temperature  at  the  twelve 
positions  in  the  tank  are  shown  in  figurec  21a-l.  The  outside 
radiation  effects  in  the  afternoon  can  again  be  clearly  seen. 

The  greatest  effect  of  this  is  at  the  surface,  smoothing  out  in 
the  pond  middle,  and  being  felt  again,  but  to  a  lesser  degree, 
at  the  pond  bottom.  This  indicates  that  absorption  occurs  at 
the  top  and  at  the  bottom,  the  greater  amount  being  absorbed 
at  the  top.  As  in  the  heating  phase,  the  time  constant  is 
smallest  at  the  pond  surface.  Figure  22  shows  the  average  tern- 
perature  variation  with  time  during  the  cooling  phase. 

5.1.3  £y^l  i.n£ 

5. 1.3.1  Equal  intervals 

Figures  23a-c  show  the  temperature  variation  with  time  of 
three  of  the  twelve  thermocouple  positions  during  the  cycling 
experiment.  Figure  23a  at  the  pond  surface,  shows  a  cyclic 
wave  with  a  10  hour  period  and  large  amplitude.  This  represents 
the  strong  effects  of  the  radiation  at  the  surface.  As  one 
moves  further  into  the  pond  the  cyclic  variation  decays  until 
it  is  no  longer  seen  at  all.  At  the  bottom  of  the  NCZ  the  tem¬ 
perature  is  nearly  constant  with  time,  revealing  the  very  stable 
nature  of  the  pond  at  this  point.  Near  the  pond  bottom. 
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a  10  hour  cycle  is  again  observed,  but  with  a  magnitude  less  than 
that  at  the  top.  This  indicates  that  the  absorption  at  the  bot¬ 
tom  is  less  than  that  at  the  top. 

5. 1.3. 2  Unequa]^  Ifiter^als 

Figures  24a-l  show  the  results  of  cycling  tests  run  in 
the  form  of  5  hours  on  and  3  hours  off.  The  same  type  of  beha¬ 
viour  as  described  in  section  5. 1.3.1  appears  in  these  figures. 

The  experiments  were  run  by  heating  the  tank,  then 
performing  the  5  on/5  off  cycling  beginning  with  the  hot  tank. 
Then  the  5  on/3  off  cycling  was  begun  immediately  after.  Thus, 
one  would  expect  to  see  the  temperatures  decrease  as  the  maximum 
overall  temperature  decreases  (see  §  2.5),  then  decrease  again 
as  the  average  heat  input  is  again  decreased.  The  series  of 
data  taken  did  exhibit  this  behaviour,  but  the  temperature 
changes  were  very  small  and  occurred  over  long  periods  of  time. 
Thus,  it  is  not  possible  to  illustrate  this  behaviour  with  the 
data  shown  here. 

5 . 2  Numerical 

5.2.1  H.e£t  2n£ 

Figures  25a-c  show  comparisons  of  the  numerical  results 
with  the  experimental  data  for  three  times  in  the  heating  phase; 

5  hours,  24  hours  and  48  hours.  Although  the  shape  of  the  cur¬ 
ves  agree  in  all  cases,  the  actual  temperature  difference  between 
experimental  and  numerical  increases  with  time.  The  differences 
are  due  to  the  simplifications  in  the  numerical  model  in  descri¬ 
bing  the  absorption  of  light  in  the  pond.  Also,  it  is  expected 
that  the  experimental  model  would  reach  higher  temperatures  than 
predicted,  as  the  additional  radiation  received  through  the 
window  was  not  included  in  the  numerical  study. 
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5.2.2  Co£l  T^n q 

Figs.  26a, b,c  show  similar  results  for  the  cooling  phase.  The  same 
comments  apply  to  these  as  those  made  in  section  5.2.1.  There 
is,  as  shown  in  figure  28,  a  difference  in  slope  at  the  pond 
bottom.  From  this  one  can  conclude  that  the  bottom  losses  are 
not  correctly  modelled. 

5.2.3  £y£llf>£ 

The  results  from  the  cycling  experiments  appear 
graphically  in  Figs. 27a, b.  Note  that  the  greatest  difference 
between  the  two  models  occurs  at  the  pond  surface.  This  supports 
the  conclusion  that  the  difference  results  from  : 

I  (1)  Pond  properties  not  modelled  -  especially  the  plastic  cover. 

I  (2)  Radiation  input  from  the  windows. 

(3)  Simplified  absorption  model. 

! 

I 

5.3  Concentration  measurements 

'I 

'  Figure  23a  is  a  concentration  profile  taken  during 

'  the  filling  of  the  pond.  One  can  see  the  di stinct 1 ayers , 

'i  that  is,  the  evidence  that,  during  filling,  the  layers  did  not 

j  become  mixed. 

I 

i 

'  Figure  23b  shows  the  concentration  profile  during  the 

i 

I  tests,  verifying  the  stability  of  the  pond. 

1 

‘  Figure  28c  shows  the  concentration  profile  after 

completion  of  all  the  tests.  A  gradient  is  still  present, 
verifying  the  stability  of  the  pond,  but  there  is  clear  evidence 
of  the  upward  diffusion  of  salt.  There  was  initially  68.4  kg 
of  salt  added  to  the  pond.  Taking  the  area  under  the  curve  of 

( 

j  figure  28c  to  compute  the  amount  of  salt  yields  63.7  kg  salt 

{  remaining  in  the  pond.  This  represents  a  1%  salt  loss. 
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Some  salt  was  definitely  lost  from  the  solution  by  deposits  on 
the  tank  bottom,  and  also  through  small  leaks  in  the  pond  near 
the  thermocouples.  But  the  concentration  measurements  are 
highly  uncertain  and  thus  no  quantitative  analysis  can  be  made 
concerning  the  salt  loss. 
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6.  CONCLUSIONS 

Both  the  experimental  model  and  the  numerical  model 
exhibited  behavior  characteristic  of  solar  ponds.  That  is,  the 
stable  region  is  at  the  bottom  of  the  NCZ  with  absorption  in  the 
upper  and  lower  convecting  zones.  There  is  heat  storage  at  the 
bottom  evidenced  in  both. 

The  numerical  model  was  found  to  be  very  sensitive  to 
the  depths  of  the  three  zones,  which  are  chosen  a  priori.  It  is 
also  very  sensitive  to  the  model  used  for  absorption  in  the 
various  layers.  That  is,  to  effectively  use  the  model,  one  must 
know  how  radiation  will  be  absorbed  in  the  pond  being  modelled. 

The  experimental  model  acted  in  quite  a  "one-dimensional" 
way,  as  supported  by  the  agreement  between  the  numerical  model 
and  the  experimental  data.  To  ensure  the  same  absorption  charac¬ 
teristics  as  in  a  full-size  pond,  one  need  only  respect  the 
agreement  of  water  quality  and  construction  materials, if  actual 
sunlight  is  used. 

The  problem  becomes  more  difficult,  however,  when 
artificial  light  is  used.  The  spectrum  of  the  artificial  light 
is  not  exactly  the  same  as  the  spectrum  of  sunlight  and  thus  the 
absorption  characteristics  are  also  different. 

Finally,  from  the  modelling  attempted  in  this  work, 
one  can  conclude  that  the  modelling  of  solar  ponds  both  experi¬ 
mentally  and  numerically  is  not  only  possible  but  feasible  and 
should  be  investigated  further. 


25 


7.  RECOMMENDATIONS  FOR  FUTURE  WORK 
7 . 1  Numerical  model 

The  numerical  model  could  be  improved  by  better 
modelling  the  pond  loss  (i.e.,  temperature-dependent  heat  loss 
coefficients).  The  validity  of  the  analysis  of  the  important 
dimensionless  parameters  could  be  checked  using  the  numerical 
model  . 


It  would  be  interesting  to  compare  the  results  of  the 
model  with  full-scale  field  tests. 

7 . 2  Experimental 

A  longer  test,  with  different  cycling  times,  could  be 
run  to  observe  the  phenomena  described  by  section  2.5. 

The  experimental  data  could  be  compared  with  full- 
scale  field  tests. 

If  artificial  light  is  to  be  used,  an  in-depth  study 
of  the  absorption  characteristics  of  the  light  in  comparison 
with  sunlight  should  be  made. 

The  concentration  probe,  as  described  in  section  3.1.4 
could  be  further  improved  to  reduce  experimental  uncertainty 
and  make  it  easier  to  use. 
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number 
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11.5 
13 

14.5 
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18 
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24 
27 
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APPENDIX  1  -  PROGRAM  LISTINGS 


(a)  Numerical  -  cooling 

(b)  Numerical  -  heating 

(c)  Numerical  -  cycling 

(d)  Program  to  solve  tridiagonal  matrix 

(e)  Plot  -  temperature  vs  height 

(f)  Plot  -  temperature  vs  time 
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F.OHMA'K  IX, 'ENTER  iHClCENT  HAL'IATIUn'1 .  . 

HtAD(5, 15)HI 

'  IIHU  IS  THE  EFFECTIVE  EXTTNCTTDN'CUEKkiCIEnT  IN  1/M  - 

UHU=.  1 

Hlrt  IS  THE  REFLECTED  RADIATION 

itlP=.6*FRAC*!II  _  _ 

7IK  =  0, 

PHIAMB  IS  THE  AML'IENI  TENPERATURE  .  . 

....  PHI  ALH  =  35./TZERU .  . . 

C(.'iE=-HH*  (i:fc  PTH  +  xTUP)**2/{TC*XTUP*100,  ) 

“  Al(l)3(-0*rELT)/C2,»DFLX}  - - 

Ub  IS  THE  UCrrUN  LOSS  COEFFICIENT  IN  v«/M2K 
Ul)  =  0. 

CTV.Gs-Ub*(XTUP  +  i)F  PTH)  »»2/(TC»XbUT  ♦100.  ) 

bl(n  =  l.t(C*n.TLT)/(2.*LELX)-Cr  lF*PFLT/2.  -  - 

t)l  (.0  =  1  .-(P*DELr)/(2.*LFDX)-Cl  wU^liEI.  1/2. 

Cl  (  1  )=  I. 

cur.  )  =  (P*  PELT) /(2.»DtLX) 

AU’0  =  C. 

nc  100  i=2,i.-i 
0=1  . 

Aid  )  =  (-0  +  CKLr)/(2.^CEL/**2) 


ri4(J()  dl(I)  =  l  .  +  (L*r.!iLn/(nr  LX**2)  - -  - 

dbOO 

H60<l  100  CC'tHN’If'. 

«7oo  c  ri“t.  Lunt>; 

d7bO  11^0=0 

h8oo  d.c;  200  j  =  i,it‘'ax 

y900  Tst  +  LEI/r 

d957  ITAusllAU+l  _ 

9000  C  CALCUl.ATt;  01  ! 

9100  C  C/^LCOLAft  HjEjAfiL  ti 

9200  H  =  +  ♦‘♦2»i'rllAf.t>/(  lC»XTUP*lUO.)  +  CUEl'rHfXTUF)** 

9  300  A  2»HI/ (TC*  l2EK'n*xTi:F*100.  J 

9400  C  IKxr  (-U.M'J*  (  (LIEFTH  +  XmCiI  )  /  too.  )  )  )  ♦  I  OtPT'l  +  XTaP' 

95uo  c  I  )**2/(rc»r 

9h00  C  H  ZEKUPXIJF+IOO. } 

X _  H  =  0,  „  _  ^  _  _ _ 

9700  [iij  HOO  IK  =  l,iJ 

99  0  0  01  t  IK  )=tii  ♦t.  XF(-tMt(*(FL.CAT(  T  h-t  )«(  OtF  f  9+ aTM(M  )  /  I  0</.  1 

9900  ,M  +riIFXt  Xt'(-U.li]»  (  (  i-FLUAT  (  J  K-1  )  )  *  ( iJKFlH+Alul' ) +XrlUl  )  / 

lOOOO  ('<  100.3 

iOlOO  FlIK)=(i:iP:FfHtXTnP)**2*DI(IK)FI)F'J4^DeLX/(TC*TZE9n»inn**2) 

10200  C  .|R1TF(70,20  JF(IK) 

10300  dOO  CO.OTINOFJ 

_  10400  _  _2&  „  _FOKr^M(lX, 'Fs' ,F10.'>)  _ _  _ _ 

10509  C  '  ■  FHKj  IS  ril£  Tt/iPtRAriihE  UF  Trif  GKOUiiO  '  “  ‘  . . 

■  ItfbOO . PHiG  =  2y  .  /TZfc  dU  . - 

10700  .  <;  =  Ub»(Ui::PTH+Xi<JP)-»'»2*l  lllG/(3C*xaUT»10O.  3  . 

_  lOfcOO  _ Z  +(L'tPlH  +  X  I'GP)  **2*Hl*(DtFTH  +  XTCl')  ♦  (  1  .-tXP(-ilMU*(0fFfH  +  X  r0P)/10o 

'lOdSO  '  y  '  )/(TC  +  rZFFU*XHGr*l(iO,)  -  „ 

10900 

11100  L'l  n=l  J*0FLTFFHI(2  )  )  /  (  2  .  ♦  OF.LX  ) -*■  (  I  .  -  (  Q  ♦OFLT  )  /  (  2  .  *06  OX)  ♦PHI  ( 

U200  A  1  )  3  + (LFLT/2.  )  *(H  +  n)+CCF;fc  ♦PHI  [  1  )*0fc.biV2. 

1  1  300  "  C"  .  wHl  rE(7‘j,  27  )!j(  1 3  "  ‘  . . . . . 

1  1400  27  FOPMATC  IX, iF  10.5)  . 

11500  . 0(0  )=(  (PtUtLl  3/(2.  ♦UtLX  3  +  1  .  )  ♦PHI  ( ;J  3  -  t  (  P*DFL  f  3  /  C  2  .  ♦  OFLX  )  )  *1’  

11000  B  H  I  ( 'J-I  3  +  (DELT/2. 3  ♦  (0  +  G3+C  r.-aj*PHl  (N  3  *0&  01/2  . 

~  1  1  700  C  .MKlTfc;(7n,29)D(N)  -  _  - . 

11900  2d  FUPMAKIX, ,F10.5) 

11900  DC  300  l32,iJ-t 

12000  0(  J  3  =  (i)*0F:Lr*PHI(J  4-1  )  )/(2.  *DlLX**23  + (  1  ,  -  (  0*PE[,I )  /  (  I3FLX 

12100  ■  “  D  ♦*2)  3  ♦PHl(l)  ♦  (  (0*Rt;LT)/(2.*(DfcLX^*2  3  )  3  *PHI  (  I  - 1  )  +  (  DH;LT/2  , - 

12200  K  )*(F(i3+KI)) 

12  300  C . KiHlTt.(7!.),29)0(l3  . 

12400  29  FOR’‘'Ar(lX, 'C=' »F10.5) 

12500  300  "'C.CiJl  ISUf:  . .  .  .  ~  '  - 

12600  C  CALCULATE  SCLOTICN  AT  TIP.F.  T 

12700  DO  301 

12600  A(JJ)=Al(Jd3 

129^0  .  b(.J03  =  BUJJ)  .  '■  ■"  ""  . .  . .  .  “  . 

13000  . C(dJ3=Cl  (JJ3  . 

13100  301  CCN  riMlE  . 

13200  CALL  THl()(A,H,C,i;,M,PHI3 

13300  ■  DO  2002  r=l,M  “ 

13400  TtHPli)=Pl(HI  )*TZFRO 

13500  2002  CCNTI.'IUf: 

13550  IF  I  rrAU,fJe.4)»iU  TC  55 

l3bOO  C  TPloT  NEW  SOLUTICO  “  ‘  -  . . . — 

13700  . WRlTLCbO,  193  r  .  . 

13600  19 . f.CRMAKlX, '1=' ,F10.5)  . .  . 

13900  WPITF;(5(),20) 

"14000  - 20  "FCP«AT(  rX.dX,  'X'  ,  1  2X  ,  '  SCLUTICH  •  ,  dXV*  TEHPFRATURF:  '  3  - — - 

14100  CO  500  n=1,N 

14200  WRIIF;(50,21  )X(KJ  ,PH1(F  3  ,Tf:MP(K) 

14  300  500  C(.!iTIriiJE 

14350  itm;=o.  - - - — - :: - 

14375  . 55 . CONTlfHJg . - .  . 

14400  200  . CLiPIINUE  .  . . 

14500  21  F(JPMAr(2X,Fl0.5,«X,F  10,5,«X,F10.5) 

"14600  STOP  ■  ^  “  - 

14700  e:w 


34 


aw  ITCii.KCR;  15 


1  ]  b  :  4  1  /  ,  U  1 


I 


100 

200 

300 

1. 1.iLNSIMN  A  (5b)  ,R(50)  ,C(5'))  ,rj(50  J  ,Rhl  (  50  )  ,X(5ri)  , 
Pd'Fi.SIbN  A  1  (  IbO)  ,bl  (  ICC  )  ,C1  (100)  ,H  dbO)  ,F(  loo  ) 
I'Yl't  11 

IE  HR (50) 

40  0 
500 

1  1 

I- Li-  -I A  1  (  1  X  ,  '  ti N  fv-R  blj.'iC'FlF  OF  SEGOF'b T5  '  ) 

RFaD(5,12)N 

.  bOO 

7  JO 

12 

(•CKMATdi) 

1  <  P  !■.  1  3 

8  00 

9  00 
lOou 

13 

FORMAT dX , 'ENTFR  NUMBER  OF  TIMk  STEPS') 
RFAD(5, 12) ITMAX 

T  y  i  •  t  i  •» 

1100 

14 

FCbFATdx, 'Enter  lelta  t') 

1  2  0  0 

h  £  A  P  (  5  ,  I  5  )  1 '  F  1 .  f 

13  0  0 

15 

■  e  0  8  !>*  A  1  (  r  1  0 . 5  1 

. 1400 

T  jRt 

1500 

lb 

F  ORMAT  (  t  X  ,  '  F.i  t'ER  r;fc.LTA  X*) 

1  b  0  0 

1700 

RF.AD(5, 15)CFLX 

TYPE  30 

18  0  0 

30 

FGRMATdX, 'ENTEb  FRACTICN  fJF  RAHlAIlUN  RFFLECTEP 
TT OM' J 

FROM  BO 

1900 

_ £_ 

. 2100 

. 22CO 

2300 
2400 
2500 
2fa00 
_  2700 
‘2h00 

. 2550 

. 2875 

_ 2887 

2893 

2900 

3000 

3100 

3200'- 

. 3300 

. 34  00 

3500 

- 3600' 

37  0  0 
3800 
_  3900 

400  0 

. 4100 

. 4200 

4300 
44  0  0 
4  50  0 
4600 
4700 
4725 
4750 
4800 

_  4900 

50  0  0 
5100 
5200 

_  5300 

5400' 
-  5500 

. 5600 

5700 

5800 

5900 

6000 

_ 6100 

6200 

6300 

. b400 

6500 

6600^ 

6700 

6800 

6900 

7000 

7100 

7125 

_ 7150 

7200 
7  300 
7400 
7500 


-■ . 

f-ORMAr  (  •  '  ,  ’  KriTt.r.  lir.t  CF  SOI.LIGHT  16  hOOHS')  ,  . 

H  t  A  P  (  5 , 1  5  )  G  a 

■T«P£  13  - - -  -  “ 

FORVAl  (IX, 'k-  .ftH  HMF;  CF  PARK  li4  tICURS'l 
RtAfKS,  151CFF 
KRt  99 

"“FCRyAT  nx;'*F7n:Ef<“AMFIEKrTT.7TF  DtTTTMG'  otl'r*  ) - - - 

f<fc:AO(5, 15)1506 . -  - . . . 

— X  Yl  t  .98-  .  - . . . . 

FL’RVAldX,'  EfiftR  ambient  XERR  DUkIimG  OAf-K') 

"R£AC(5,  IblTDARK  •  - . . .  .  -  . . 

PG  3000  1=2, N 
X(I)sX(I-U+DEOX 
COM  IMIfc 

“XCi)  =  0 - - - — - 

-TYFt- -luOl  . . 

.  f.CHy.AK  IX,  'EnTEH  KEFERthCE  lEM.l'tkA  1  Uht  *  ) . 

8hAU(b,  15)i:2-t:RO 

“CALCULATE  G  AND  F  - - - - 

lYFt  22 

FCKMAi  C  IX,  'ENTER  PERTH  CF  IJCZ  IN  CR  '  ) 

REAOCS,  l‘i)X'iUP 

“TYPE  2-1  - - - — - - - - 

--  FGRRATdX  , 'ENl'ER  UKRTH  OF  NCZ  IN  C*"  )  . 

.  .  REAL(5, 15)PhPIH.  -  . . . .  . . 

XYPE  23 

“TCRyATdX,  'FN'TFR  DEPTH  GF  LCZ“IH  c;"  )  - - - 

REAUCb, 1 5) XUU  r 
n=(Dr:FTH-fXlfJR)/XTCP 

p=-(D£prH  +  x'njp)/xL'ur 

“TLIGMT=.085  8Anii*6  0./CDtT'Tn  +  XTCP)»*2 - 

ILARK=.(J85  8+ftFF^6<),/(PEFTH-fXT(iP)**2-l>tLT  . . 

-CALCULA  IE  IlilT-lAL  SGLUIIGN . . 

TYRE  25 

“  FCRyATdX,  'ENTER  INITIAL  CONDITIONS  TN”  NON-DIM  FGK«  •  ) 
(VO  700  lR=l,ri 
READ(5, 15)PH1 (IK) 

CUNTINNE 

“  DG^  20C0~I=l - - - 

IFyPd)=RHld  JRTZEHU  .  . 

.Xd)=f  I.UATd-l)+D£LX  .  .  . 

CUNT  I  ONE 

“PRINT  IMITTAL  SOLUTION - - - - 

WRIT£(50,20) 

DC  600  11=  1  ,-1 

ftRllE(50,2l)X(II),RHI(lI),TEMP(Il) 

“CCNTI.-iUE  ■'  -  - - 

IMAXaFI/OAK  ITMAX)-*1)FLT  . . 

-  CALCULAXE  A(I),Bd),C(I)  .  . 

TC  IS  I’OE  TBEKclAL  CONDUCTIVITY  OF  IHE  RONn  IN  w/MK 

■“TC  =  0.596  - - -  -  "  '  -  -  “■ 

HR  IS  THE  TUP  LOSS  COEFFICIENT  IN  «/«2K 
HI  IS  THE  RADIATION  IN  W/“2 
TYPE  31 

“FOR  'ATdX, 'ENTER  INCIDENT  PADTATIO  J*  )  “  “ 

KhoPIb,  15)r)l'ip 
Hlsuli.C 
TA.'-'Rs'lSOii 

UMU  IS  THE  EFFECTIVE  EXTINCTICM  CUtFFlClFNT  TN  l/‘' 

iJRi.s.  1 

Hlt<  IS  Tile.  REEl.FCTED  RADIATIGN 
‘<F  =  o  . 


7800 
7000 
HOOD 
^  blOO 
8500 
8600 
8700 
8H00 
8000 
9000 
9100 

_  9125 

0  150 
9175 
9200 
„  9300 

0  4  0  0 

. 950  0 

9600 
_  9650 

"  ‘9700  ‘ 
9800 
9850 

_ 990  0 

1 0  0  0  0  ■ 
to  100 

10225  .. 
10237 
^10  243” 
10250 
10  2  51 

_ 10252 

loiOLI 

10400 

.  1050O  . 

_ 10600 

10  700 
10800 
10900 

_ llOOO 

moo 

. 11200 

11300 

_ 11400 

11500 

11600 

11700 

_ 11800  _ 

11900' 

12000 

12100 

_ 12200 

12300 
12400 
12500 
„  12600 
12700 

. 12800 

129O0 
13C00 
13100  “ 
13200 
13300 

_ ^13400 

i350ff 

. 1  3600 

13700 


CL  II  =-i-K*(i,'Fprn  +  XTnt>^*42/CTC*XT.lH*  lilO.  ) 

A1  (.  1  )  =  i-<j*i;kL  fj/C  2.*l;l,I,X  ) 

C  Oil  ly  I'.IL  I'.LIIUM  LU5S  CttFUC  Ihl.  f  1  .v  «/M2K 

tiLf  =  0, 

Cl(l)=o. 

C  1  (0)  =  (P*UtLTJ/(2.'H)KLX) 

Al( 

liC  100  1  =  2, W-l 
U=l. 

Al  il  )  =  (-U*0tI,T)/(2.*l;F:LX*  +  2) 

ai  11  >=l  .  +  (U*ljf:6l  j/(LlFLX*»2j 

cT.ai=-  (i!*(  xioa  +  nf.PTii j ♦♦2/(  rL*xani  ♦  1  oo. } 

oi  t  u  =  i.  +  (o*L'KLT)/(2.+i:tLX)-tcTif:*a<.:i.T//:. 

a  1  (9  )  =  i  .-iP4rie;LT)/(2.*ctiiX)-L'i..u*uhL  1/2. 

CUl  J=A1(1) 
loo  CGNTIfilJfc: 

i  iatsO.  -  ■  ^  - .  ■  ' 

. 'i'so. . 

C  . Tlf't'  i.ciuP: 

1IaJ  =  () 

uu  200  J=17ITMAX  ”  '  - 

T  =  T  +  L’t;L,T 
ITAUri  TAiJ  +  l 

_  f  IML  =  T19t:  +  CFLT 

iF(Ti‘'t;.Lr.TLiGaT)Gn  tc  loi  '  -  - -  - 

.  . IF (  I  IrtF.GE, ( fLIGH  f  +  TDAKK ) )TI«F=0.  . 

112  =  1,125 .  . 

91=0. 

lF(TlHf:.E'(i70.)i'I  =  HI«F - - - - - - 

lAfl8  =  l()Ai<K 

C  PhlAfip  JS  THE  AMBIEM'J  IfcMPERAlUKt 

_ t’HlA)Xa= rAMe/iztkt 

C  ■  CALCULAlt  Xr'“  - - - 

C  . CALCOLATE  a,F,AOI;  Q . . 

nl Ps  ,  o *Ek/»C  ♦  H I . 

_  101  H  =  nR  +  (iiFPTH  +  XI'UP  )  *♦  2*FH  I A  96/ (  Tf*X  1  LiP*  1  00  ,  J  +  (  nEETH  tX  TOP  )  ♦  ♦ 

'  '  A  “2*riI/(TC*T2£KO*XTCP*lO0.)  - —  -  ^  _  _  — 

G  +illP*( 1 ,-.o4EXP(-0H0*( (DEPTH+XHUT-6. )/lOO, ) ) ) ♦ COEPTH+XTUP 
1  )»«2/(J'C*T 

_  ___H_ZE(<0*XrijP*109,  ) 

DG'  800  'IKsl  ,0  -  - - - - 

. UK  TF  )s:Ht»EXP(-IMtJ*(FLLATt  IK«1  )*lL»tPlh  +  ATUP)  )/100.  J 

. M  •HtlH*,t>*EXP.(-UM04  (  (l-f  LGAKIK-l)  )  ♦  tC  tP  l  H  +  XTUP  ) +XBUT-6  ,  )  /  . 

_ f)!  100.) 

F  (  lK)s:Tr)E'PTHfXTOF)**2*DI(IK)4m-HJ4DFI7X/CTC»TZERO*lOa**2r 
C  WKlTE(7o,26)F(lK) 

800  CCNTIMUE 

_ 26 _ t.GHNA  KlX.  'E='  ,F10.5) 

C  PhIG  1.5  THE  TEMPtHATURE' OF  THr  GHUUivC - - - 

. PHiG=  I5vi/TZERU  . 

G=ll64CUEPTil+XlUP)»*2*PHIG/(TC*XbOT4  100.  )  +  (DF.PTH  +  XTOP)*» 


1  3800 
1  3900 
14000 
14100 
_d4200 
14300 

29 

—  300" 
C 

14  4^/0* 
14500 

101 

14600 

14  700 

14800 

14900 

2002 

14950 

C 

.-FXF(-OMO*(UEP7H4-XTQP)/lt)0,  )7/(TC*T2EFO*XBUTMOO.  ) 
■+(DtFT’H4-XTljP3*<‘24'.6»FHlAhB*Hir»EXPl-UMU»(XB0T-6.1/100.  )/( 
TC*TZt80*XBUT *100. ) 

Ul  1  J=(Gl*UELr*PMI(2)  J/(2.*DtLX)  +  (  1  ,-(0*DfcLT}/(2,*0ELX)4PHl( 
nj  +  (PELT/2.  )♦(;(♦  ID+CrjkFKPHI  (1  )*UELT/2. 

~XT<TTE(70,27)D(1)"  - - - - - — - 

■  FliHrfAIdXy  'Uls' ,F10.5J  . 

-a(H)  =  (  (P4'|jELT)/(2.*DEI.X)  +  1.  )  46111  (Id- ((P*DELTJ/12.*6ELX)14P 
MI  (W-1  >  KMELT/2.  >*(G  +  G)+CTWU»PII1  (N7  *U£LT/?, 

“rtPri T E  ( TO-,  2  8  )  D  CM  ) -  - 

FOkMAK 1 X, 'MNs* ,F10.5) 

DC  300  1=2, N-1 

_D(I)  =  (U*rjFLT4PMr  (I  +  17)/(2,*DEI.X**2>4ll.-(O»0ELT)/(DELX 

♦  ♦27J*PHraiT(  (u4CEn,T)/(2.4  (CELX»»n  )  J*PHI(r-l  )  +  (DELT/2, - 

-O  +  CF  (1)+E(1J  ) .  . . 

«lkITEC70,29)D(l)  . . . 

FCK^!A1(  IX,  «D=' ,F10,5) 

CCHnoUE -  - - 

CALCULATE  SOLOl ICM  AT  TIME  T 

DC  301  JJ=1,N 

A(JJ}=A1(JJ} 

'B  (  J  o  )  s  til  (  JJ  )" — - — - 

C(Jd)=GHJvJJ  . . . 


15000 
1510) 
152CO 
1 5  i  0  0 
15409 
1550  0 
15600 
15700 


_  CALL  TRIO(A,H,C,D,N,PHn 

no  2  on  2  “I = r/H - - 

TEMP(1)=PH1 (I)4TZEH0 
CUM  I  HUE 

IF(  I  fAU.EE.DGC  70  55 

r^RPIT  MEW  SCLUTICfl  ^  - - 

-iRllfc  (50, 19)  I 

F.Gk  ik  t(  lA,  'is'  ,F1Q,5) 

«RlTE(5O,20) 

I.CRmak  tx,8X,  'X'  ,12X,  'SCLIITICN'  ,8X,  '  TEMPER  ATUPE  '  ) 
CU  500  6=1,6 

«Rnt(50,21  )X(k)  ,PHl(K),TtVP(k) 

CDi'Ti  ii  OE 


ib750 

1577b  C  55 
15«00  i^OO 

15900  21 

16000 
IblOO 


lTA!i  =  0  -  oo  - 

Cv  M  T  1  r-j  U  h, 

CCfJl  l.tijf. 

FuH  M A  1  I 2X , F 1 0 . 5 , 6X , FI C , 5 , OX , f 1 0 . 5 ) 
i)  T  C  P 

F  O’  1> 


10('  Sijt.-hbi' 1 1  t.h  I'Kli:  ( ,Li  ,C  ,r  , ,  i- ia  ) 

h  i; '}  C 

3  3  GO  C 

iib't  Olr  fc-SHKl  <>(  I  )  ,r;(  u  ,<-■(  1)  ,t  n  )  ,PH  H  1) 

J2(in  C  uGLvt  FON  TOF  .‘Oil.ijr  1(.::  VtCTd!-  n\  uCol-LL  G-sF.FF  ALiCGOHl in’' 

Jj'iij  C  t  .lKciT  o  .KI'.t  “  f. L 1  '  1  AT F  L  JV’F'F  i  lA>iC'«Alj 

i’iC'i  LC  1  1  =  2, 

3b')-'  F  ( I  )=t  (1  )-(■  (  1  )  »»(  ?'1  )/l-  (  1-1  ) 

3n-J(i  r  =  n-FC  L)  IZ-rd-l  ) 

3700  <:(U  =  n. 

3>lOO  1  C(  -Frii. 

3y')0  C  GFtU:;!  ')Sr>EFl>  -  r.f'l'FF '-li'h  .'il.'Ll'TTUfl  VkCTuF-  • 

40v;'J  pt:i  (-i)=i-(’  J 

41'jO  I'f.  2  .l  =  i,ii-l 

4200  0=0-1 

4300  2  FifI(J)  =  (D(o}-A(o)»!'Fl(J-(  !))/»(. J) 

440'‘  C  SCLoTKm,  CDi'.fLtl'F 

f  4boO 


(e)  ■  ■ 

P  'j  u  R  H  M  5  R  0  y  N 
L  0  j  I  V  i_  +  i  Y  ti  'i 

Y  :  PEYY  !0i  I  YEC'-  1  >.  J  VEC(  4  >,  TC  1  r  >  ,>((:  17  >,  TH(  15  ) 

YsTh  /3  .  .  e..  .  R.  .  1  1  .5,  13  .  -  14  .  5.  15  ,  ,  18,  ,  21  .  ,  24.  .27  .  .  30 
33  .35  . 3?  .  . 42  .  45  ,  / 


EPiLCULhTT  ONS  WITH  OUTFIT  LIST  [1]'/ 

■  PLOTTER  C23  ' / / 

C  ■:  E  ~  '  :  0  1  f-  ,  £  0  p  T 
FORrisT<- :  5 

:  F  :  i  0  F  T  £  Q .  1  :<  c  ft  l  l  assign  'S.’bsown.out';- 
!  F:.  !  OFT  .  EG  ,  n  00  TO  5 
TYPE  1012 

fO.RN^s  r'  -'  TOTftL  NUMBER  OF  DPiTh  POINTS  ->  '  $■  > 

H  C  0  c  P  T  1  0  I  0  .  N  P  T  0  T 
nOLE  =  ^' 

T  V  p  p  7  ~  7 

FC'RMkT  ■' D  0  YOU  UftNT  STARS?!  Y/N  3'  > 

I  F':  YE3v  J  3MODE  =  7 
£  D  £  S  - 
I  P  L  0  T  =  F; 
t1  Ii1ft=0 

I  F£  I  0=  T . EQ,  1  )  GO  TO  5 
I  V  ?  E  1  0  1  1 

FOiTtlftT''-'  HURSER  OF  POINTS  TO  SKIP  -> 

ACCEPT  1  0  I  ■■> ..  H  S  K  .1  F 
IF(  NSl'.  IF  .  EQ  .')  >  MSKIP=1 
f!CUP.VE  =  HPTOT  ^HSK  IP 
NPLOT^O 
I  1=1 

DO  25  1=11,4 

CALL  GETVAL  C 0  ,  I V£ C , 1 S T ft T > 

I  F<:  I  ST  AT  .  LT  .  0  )  GO  TO  099 
1 FC  IDEE. NE . 0  )  GO  TO  21 
TYPE  15 

FORMATC/-'  T-REFERENCE  ->  '$) 

IDEE=1 

ACCEPT  20,TREF 
FORflAT’:  F  1  0  .  A  ) 

J  V  E  C  <  I  >  =  I  V  E  C  <  1  > 

C  0  (I  T  I H  U  E 

I  FC  I  OPT  .  EC!  1  )  TYPE  9000,. JVEC 
FORMAT''.  1  :>!  ,  4  1 1  0  > 

A  VE  =  C  VECC  i  )+.j  VECC  2  >  +  .JVEC(  3  )  + JVECC  4  )  D-M 
C-  0  3  0  1  =  1,1  7 

CALL  GETVAL  C  C- ,  I  VEC  ,  I  STAT  ) 

TCI  )  =  T  R  E  F  + 

1  C-3  .0.3  7SaRTCb4.4S  +  0.02564*CFL0ATi:  I  VECC  1  )>)>>/0. 01282 

CONTINUE 

I PLOT= lPLOT+1 
I FC  I  OPT . EG . 2  )  GO  TO  2000 
IFCMIMh  NE.O)  go  TO  2001 
E  M  A  =  ; 

.cm N  =  "''  :  ;• 

, :  N  p M I N 
TO  2 0 2  1=1.17 

IF  .I  .E  ,!  1  4  OR  .  I  .  EQ  .  16  >  GO  TO  2002 
,:m  I  N  =  ftm  H  i;  :!:i  i  n,  Tc  I  > ) 
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n  r;  ^  M  =  (< !'!  X  1 '  X  M  A  x .  T  ■:  I  >} 


2  0  c  2 

:OXT  l;r.i£ 

4  P.  I  T  t  •  .  2  .  3  5 

0  ? 

Z  5 

F  0  P  n  A  T  ■:  I  ; 

‘0  F  A 

- 1  / 

'v  RITE  '■  2  j  4  0 

\ 

4 

■  0  p  ri  Fi  T ;  1 H  ;■ 

2  V  0 

lOYT  !:tU£ 

if; I  CP  t , eg  - 

1 

G  0  T  0  3  0  0  0 

IF'  NSt:  IP  .  ES 

.  i  .* 

GO  TO  2003 

•  F  i  ;i  c  0  ■  ;  I  p  L  a 

T  ,  H 

SK  I'^  )  ME  .  1  ;■ 

GO 

T  0  '5  V  '■J  0 

2  0  0  3 

I  T  =  0 

0  0  2  0  1  =  1 

*17 

I  F>  I  .  EO  ,  1  4  . 

I  T  =  I  T  1 

T  H ;  I  T  ;  =  T  ;  I 

OR  . 

I . E3 . IS  >  GO 

TO 

2  0  1  ii* 

2  0  10 

C  0  N  T  I  li  ij  E 
CALL  LILI  ( 

15  .. 

TH , K , 1 , RODE 

.  HC 

URVE, 3  ,  1  ,  1 

r. 

I  P  L  0  T  =  0 
i)PLOT=H?LOT 

+  1 

C 

3  0  0  0 

DO  50  I = i . 2 

7' 

call  GETVAL 

C  0 

..  IVEC,  ISTAT 

) 

I  F  (  I  0  ?  T  .  £  0  . 

I  > 

GO  TO  3001 

I  F  C  N  P  L  0  T  .  E  Q. 

N CURVE)  GO  TO 

52 

3  0  0  1 

I  F  I  S  T  A  T  .  L  T 

.  0  ) 

GO  TO  999 

5  0 

C  0  H  T  I  ri  0  £ 

S  1 

CALL  GETVAL 

C  0 

, IVEC, ISTAT 

) 

I  F  (  I  S  T  A  T  .  L  T 

,  0  > 

GO  TO  999 

I  FC  I  AS  SC  I VECC  1 

) > . GT . 1 00  0  ) 

TO  51 

J  VECC  •  >= I V£ 
11  =  2 

C(  1 

) 

GO  TO  10 

5  2 

IFC  ISTAT  ,  LT 

,  0  > 

GO  TO  999 

CALL  GETVAL 
GO  TO  52 

C  0 

, IVEC, ISTAT 

> 

9  99 

I oe8=o 

TYPE  55 

5  5 

FOPdATC’  HE 

RUM  [Y/N3 

; 

'  $  ) 

I F( YESC ) j  UO  TO  5 
I  F<:  I  OPT  .  HE  .  1  >  GO  TO  999  9 
i.!RI  TE'  2.  35> 

CLOSE  CUNtT=2,D13?*'PRlNT' > 

STOP 

END 


9999 
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-'HS 

'  I  iEr; ; ; : :i  ;  vsc  :  i  j  vec(  4  t^:  1 7  • .  1000  > .  tni  1 000  > 


'  TGT^L  OF  Oii'N  POINTS  -> 


:  J=  22773  «^FL0PT.  i-i  > 


0  0  2  2.  :  -  M  4 

CA_L  •F'-JiiL  0  .  I  VEC  ,  1  STF.T  > 

:  17  7;-'  L  '  O  ;  GO  TO  90? 

;c';:7  ■■(£  o;  go  to  21 

FCFfV--.,'  T-FE'EFEttCE  ->  '$> 


0  •: :  E  p 


:•>  .  -"REF 


PROBE  MUHBER  .ENTER  tOO  POP.  AVERAGE, ->  '»> 


-■."EP"  2  •'■2  HP  POSE 


V  EG  !  ,  =  I  V  £  C  (  i  ‘1 


AVE  =  '.!VEC<l>+jVEC(2:+jyEC(7>  +  JVEC':4>;'/4 


:•  0  3  0  :  - 1 . 1 7 

CALL  GETVAL  T  <>  ,  I  V£  C  ,  I  S  T  A  T  > 

T'.  r  >  =  *P£F  + 

1-3  '.;7  +  SaRT-;  34  .  43+0 . 02564*':  FLOAT!  I  VEC:  1  >  >  >  )  )/0  .  VI  282 
COifT  Iti'JE 

I T« I T+  1 

!  F!  NPROBE  EQ.  1 00 >G0  TO  500 
TN:  I  T  >  =  T( NPROSE  ) 

GC  TO  3000 
S  U  !1  = 

0  0  7  0  J  = 1  .  1  7 
sijn  =  suM+T<  J  ;< 

COST  I  HUE 

SUA  =  ?L)ri-T<  1  4  )-TC  16) 

TN':  I :  .'-sun.'  15 


[  I 


DO  5 C  1=1.27 

CALL  GETVAL  ( 0  ,  I VE C , I  ST  A T ) 

Iff  ISTAT  LT  .  0  )  CO  TO  999 
CO;iT  I  HUE 

CALL  GETVAL  < 0 , I VEC . I  ST  AT  ) 

IF(  ISTAT  LT  .0  )  GO  TO  999 

IF!.  lACS!  I  VEC<  1  )>  .GT  .  1000  )  GO  TO  52 

J VEC!  1  )= I VEC<  I  ) 

I  1  =  2 
GO  TC 


!F  ;  ;T-T  LT  V)  GO  TO  999 
CALL  GET'AL  f  0  ,  I  V£  C  ..  I  ST  A  T  ) 
GO  TO  32 

t  r  .  7  r»  -  ^ 
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APPENDIX  2  -  APPARATUS  INFORMATION 


Insulation 
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PITTSBURGH 
CORNING 
EUROPE  SlA. 


[p;?;]OP3=5eiET 

Deux  types: 

1  -  FOAIVIGLAS®  T2  COMPOSITION  : 

principalement  destine  aux 
applications  pour  lesquelles  la 
conductivite  thermique  repre¬ 
sents  le  facteur  predominant,  temperature  D’EMPLOI  : 


=s  PHYS301LJE3 


Verre  cellulaire  alumino-silicat6 
d’une  composition  specialement 
^tudi^e,  totalement  inorganique,  sans 
addition  de  liants 

-260  a  -I-  430*C 


POINT  DE  RAMOLLISSEMENT : 

2-F0AMGLAS®S3. 

applique  plus  particulierement  absorption  d’EAU  : 
lorsque  les  exigences  en  (A  I’immersion) 

mati^re  de  resistance  e  la 
compression  sont  severes.  hygroscopicite  : 

PERMEABILITE 

ALAVAPEURD’EAU: 

Deux  isolants  FOAMGLAS® 

pourront  done  etre  specifies 

pour  les  applications  dans 

les  domaines  du  batiment,  Qg  vapeur  D’EAU  : 

de  I'industrie  et  du  froid. 


Point  de  ramollissement  du  verre 
environ  730“ C 

Nulle  a  I’exception  de  la  retention 
momentan6e  de  surface 


Nulle,  selon  norme  AFNOR-NFX41.001 
(Essai  du  C.S.T.B.  -  Paris) 


//  Forschungsinstitut  fiir 
Warmeschutz,  Allemagne 
(Dr.  Cammerer) 


CAPILLARITE : 


RESISTANCE  AUX  ACIOES : 


RAsiste  A  tous  les  acides 
communAment  employAs  ainsi  qu'A 
leurs  vapeurs 


INCOMBUSTIBIUTE : 


ClassA  M.O. 


STABILITE  DIMEMSIONNELLE : 


Parfaite 


ISOLEMENT  ACOUSTIQUE 
BRUT  AUX  FREQUENCES 
AUDIBLES : 


28  d3  pour  une  Apaisseur  de  10  cm 


Dimensions  en  mm. 

Dimensions  en  mm. 

Longueur 

Largeur 

Epaisseur 

Longueur 

Largeur 

Epaisseur 

300 

450 

40, 45,  50,  60. 70.  80,  90, 100. 120, 130 

300 

450 

40,  50,  60, 100 

600 

450 

50,  60. 70,  80.  90. 100. 120. 130 

600 

450 

50,  60, 100 

marquage  rouge 

marquage  bleu 

M 

SI 

M 

SI 

MASSE  VOLUMIQUE 
(tolerance  ±  10%) 

'  125  kg/m’ 

125  kg/m’ 

135  kg/m’ 

135  kg/m’ 

CONOUCTIVITE 
THERMIQUE 
(tolerance  ±5%) 


*-20"C:0,033l(cal/m.h.*C  i  - 20“ C:  0,038  W/mK 
0“C: 0,036 l(cal/m.h.“C  0“C:  0,042W/niK 

+20“C:0,039lical/in.h,“C  +  20“ C :  0,045  W/m K 


»-20"C:0,035l(cal/m.h.“C  i  -20“C:  0,041  W/mK 
0“C: 0,038 kcal/m.h.“C  0“C:  0,044  W/mK 
♦  20“  C :  0,0  4 1  kcal/m.h.“  C  ♦  20  “  C :  0,048  W/m  K 


•RESISTANCE  A 
lA  COMPRESSION 
(valeur  moyenne) 

'  5  kg/cm’ 

490  kPa 

7  kg/cm* 

680  kPa 

RESISTANCE 

A  U  FLEXION 

4,5  kg/cm’ 

440  kPa 

5,3  kg/cm’ 

520  kPa 

MODULE 
ffElASTICITE 
(A  U  FLEXION) 

10.000  kg/cm’ 

980  MPa 

12.000  kg/cm’ 

1.180  MPa 

COEFFICIENT 

DE  OIUTATION 
LINEAIRE 

1 

1  8,5x10V“C 

8,5x10  '/K 

O.OxIO’/'C 

8,5x10  ’/K 

CHALEUR 

SPECIFIQUE 

1 

;  0.20  kcal/kg*C 

0.84  kJ/kg  K 

0,20  kcal/kg’C 

9.84  kJ/kg  K 

DIFFUSIVITE 
THERMIQUE  A  0‘C 

4.0  X 10"’  cm’/sec. 

4.0  x  10"’  m’/sec. 

4,0  X 10  ’  cm’/sec. 

4,0  X 10'’  m’/sec. 

*Le  coetticient  de  s6curH6  sera  choisi  par  le  bureau  d'itudes  en 
torjction  de  rappUcation.  Une  valeur  de  3  est  triquemmeni 
adoptie.  Saut  indications  contraires,  les  propri6t6s  du  verre 
celiuiaire  F0AMGLAS9  sent  donnies  i  ia  temperature  ambiante 

FABRIQUE  PAR 

PHTTSSURGH 
CORNING 
EUROPE  SA 


PITTSBURGH  CORNING  EUROPE  S  A 
Avenue  de  Tervuren  36  -  Boite  19 
B-1040  BRUXELLES  Belgique.  T6I  (02)  735  90  36 
Telex  22277  PiClO  B. 


et  correspondent  aux  methodes  d'essais  ASTM  N°  C303.  C240, 
C 165.  C  177,  C  203.  C  35S  et  E  136. 

Pour  Visoiation  des  londs  de  reservoirs,  on  se  reierera  i  ia 
specitication  i  22  -  edition  1978. 


I 

e 

J 


f 


t^fOAUSliS  u  K  tMt  /«  0ipBS*n  ftff  ftitt  I’/iii  tt  im  ftuun 


FIGURE  1 

RADIATION  PHENOMENA  IN  A  SOLAR  POND 
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FIGURE  2 

LIGHT  ABSORPTION  IN  WATER 
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FIGURES  ;  RADIATION  nELD  ISOLINES 
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TEMPERATURE  MEASUREMENT  CHAIN 
FIGURE  7 


Volts 


g/mi 


POND  FILLING 


RADIATION  EFFECT  ON  TEMPERATURE  MEASUREMENTS 

FIGURE  14 
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EXPERIMENTAL  DATA  -  HEATING  -  HT.=30.CM 


Cj  2338030)  3aniV83dW3i 


EXPERIMENTAL  DATA  -  HEATING 


EXPERIMENTAL  DATA  -  HEATING  -  HT.=3.CM 


FIGURE  17 

AVERAGE  Temperature  during  heating 


RGURE  19 

OVERALL  LOSS  COEFFIENT  DURING  HEATING 
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EXPERIMENTAL  DATA  -  COOLING  -  HT.=27.CM 


FK3URE  21c 
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EXPERIMENTAL  DATA  -  COOLING  -  1 3. CM 
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EXPERIMENTAL  DATA  -  COOLING  -  HT.=9.CM 


EXPERIMENTAL  DATA  -  COOLING  -  HT.-3.CM 


FIGURE  22 
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AD-«t3e  703  THE  $ALT-M*01ENT  SOlU  fONOlUI  VON  MMUN  INST  FO* 

FLUID  OTNANICS  HHOOC-MINT-aCNESE  KELOIUNI  S  T  RROVN 
FED  83  E0A80IRS3-4  4FOS«-t2-0>OI 

I/NCLASSIFIEO  F/0  10/3 


NL 


r 


MICROCOPY  RCSOLUTfON  TEST  CHART 

NATION/M.  tUftCM/Of  sr*NO*M>S-'l»C9*A 


/ 


8CSii;«!(98888*'‘ 
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FIGURE  24  j 


POND80.  PROBE  3  CHT.  =9.CM> 


POND80,  PROBE  1  <HT.=3.CM) 


FIGURE  25  a 


FIGURE  25c 


FIGURE  26  b 


CYCLING  -  BOTTOM  OF  NCZ  CHT.=t2.CM) 


FIGURE  28  b 


bottom 


FIGURE  28c 

concentration  profile 


CONCENTRATION  (g/ml)  ®  P^obe  measuremenis 


